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Abstract 

Using Monte Carlo simulations we study a coarse-grained model of a water layer confined in 
a fixed disordered matrix of hydrophobic nanoparticles at different particle concentrations c. For 
c = we find a first-order liquid- liquid phase transition (LLPT). For c > our simulations are 
consistent with a LLPT line ending in two critical points at low and high pressure P. For c = 25% 
at high P and low temperature T we find a dramatic decrease of compressibility Kt, thermal 
expansion coefficient ap, and specific heat Cp. Surprisingly, the effect is present also for c as 
low as 2.4%. We conclude that even a small presence of nanoscopic hydrophobes can drastically 
suppress therodynamic fluctuations, making the detection of the LLPT more difficult. 

PACS numbers: 64.70.Ja, 65.20.-w, 66.10. C- 
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Many recent experiments investigate the behavior of water in confinedgeometries [1] for 



its relevance to nanotechnology, e.g., filtering water in carbon nanotubes 



2j, and biophysics, 



e.g., intracellular water {3]. An interesting property of nanoconfined water is that it remains 
liquid at temperatures where bulk water freezes. The present technology allows us to ob- 
serve bulk water in its liquid phase below 0°C if quenched very rapidly (supercooled), but 
ice formation cannot be avoided below T# = — 41°C (at 1 atm). Interestingly, a number of 
theories and models predict a peculiar thermodynamic behavior for bulk water below Tjj, 
with a liquid-liquid phase transition (LLPT) ^-6|. Although studying nanoconfined water 
could shed light on the phase diagram of deeply supercooled water, experiments and simu- 



lations {2I show that fluid-fluid phase transitions in a confined space can differ from those in 



bulk water. Several studies using specific geometries, e.g., slits |8-ll0j] or disordered matrices 



of disks or spheres 
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12], have clarified some aspects ^j u ^ l eave open questions about the 



thermodynamics of supercooled confined water |l|, |9_|, |l3j, |l4j . 

It has been proposed that supercooled water forms highly structured regions in the hy- 
dration shell of nonpolar solutes 15] , where the hydrogen bond (HB) network is weakened 



only when the size of the hydrophobic particles is above a characteristic value [16[, calcu- 
lated using free energy analysis to be ~ 1 nm Muller explained experimental results by 
assuming enthalpic strengthening of the hydration HBs with a simultaneous entropy increase 



in the hydration shell 



Here, motivated by several experiments on water in a strong hydrophobic confinement 



19j . we consider a water monolayer of thickness h < 1 nm in a volume V partitioned 



into JV cells of a square section of size \J"V jjYh. Each cell is occupied by either a water 
molecule or a hydrophobic particle. Particles can occupy more than one cell, depending on 
their size, are spherical and approximated by the set of cells with more than 50% of their 
volume inaccessible to water. Particles are randomly distributed and form a fixed matrix 
that mimicks a porous system or a rough atomic interface. N < JV is the total number of 
cells occupied by water molecules and V < V is their total volume. The Hamiltonian for 

n 

water-water interaction is |5[ 



EE U ( r ij) ~ JiV HB - Ja n > Yl ^ 

i (k,i)i 



(1) 



Here is the distance between water molecules i and j, U(r) = 00 for r < tq = 2.9 A, the 
water van der Waals diameter, U(r) = £ w [(r /r) 12 — (r /r) 6 ] for r > r with e w = 5.8 kJ/mol, 



the van der Waals attraction energy, and U(r) = for r > r c = y/VJh/4, the cut-off 
distance. 

The second term of Eq. ([T]) describes the directional HB interaction, with J = 2.9 kJ/mol, 
and the total number of HBs Nhb = ^iij) n i n j^<Jij,<jjii where rii = 1 for a water molecule 
when Nvo/V > 0.5 (liquid density, with vq = hr^) and rii = for a hydrophobic particle. A 
HB breaks when the OH — O distance exceeds r max — toh = 3.14A, because riirij = when 
the 0-0 distance r > r max = r ^2 = 4.10A (r OH = 0.96A). It also breaks if 00H > 30°. 
Therefore, only 1/6 of the orientation range [0,360°] in the OH O plane is associated with 
a bonded state. By allowing q = 6 possible states for each index er^-, we account for the 
entropy loss associated with the formation of a HB because, by definition, S aijj(Tji = 1 if 
<7ij = <jji, S aijj(Tji = otherwise. The notation denotes that the sum is performed over 
nearest-neighbors (n.n.) water molecules % and j, so that each water molecule can form up 
to four HBs. 

HB formation increases the volume per molecule, because it leads to an open network of 
molecules with reduced n.n. due to close molecular packing. We incorporate this effect by 
an enthalpy increase -Puhb for each HB, where vhb/^o = 0.5 is the average density increase 
between high density ices VI and VIII and low density (tetrahedral) ice Ih. 

The third term of Eq. (TjQ) accounts for the HB cooperativity, with J a = 0.29 kJ/mol, 
where (k, £), indicates each of the six different pairs of the four bond-indices of a molecule 
i. It gives rise to the 0-0-0 correlation, locally driving the molecules toward an ordered 
configuration fl^ . 

The water-nanoparticle interaction is purely repulsive, U wn (r) = e/j[(ro/r) 12 ], with = 
e w y/0A = 1.8 kJ/mol jl2], where r < r c is the distance between the water cell and each of 
the cells occupied by the nanoparticle. The restructuring effect of hydrophobic particles on 
water is incorporated b y r eplacing J and J a in the hydration shell with J h = 1.30 J and 

= 1.30J CT , following 20[. Because bonding indices facing the nanoparticle cannot form 
HBs, at intermediate T they have a number of accessible states larger than those facing 
water molecules, inducing an increase of hydration entropy jisj ]. 

We perform Monte Carlo (MC) simulations for constant pressure P, T, and N, with 
variable water volume V = Vq + Nhb^hb, where Vq > Nvq is a stochastic continuous 
variable that fluctuates following the MC acceptance rule 21[. We simulate systems with 
jV < 1.6 x 10 5 within a fixed matrix of spherical nanoparticles of radius K— 1.6 nm, with 
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nanoparticle concentration c = (jV — N~)/jV = 2.4% and 25%. We repeat the analysis for 
R = 0.4 nm. For c = 0, the model has a phase diagram with a first-order LLPT, between a 
low density liquid and a high density liquid, starting at P ~ 0.2 GPa for T — > and ending 
in a critical point at T ~ 174 K and P ~ 0.13 GPa Q. 

We find that for c > the liquid-gas spinodal is shifted to lower T and the line of 
temperature of maximum density (TMD) is shifted to lower T at low P and to higher T at 
high P, with respect to the c = case, reminiscent of results for other models of confined 
water fiol . \\\ . We find stronger changes for increasing c (Fig. [T|). 

Further, we next find that confinement drastically reduces volume and entropy fluctu- 
ations at low T. To quantify this reduction, we calculate volume fluctuations, entropy 
fluctuations, and cross-fluctuations of volume and entropy, and analyze the associated mea- 
surable response function, respectively, isothermal compressibility Kt, isobaric specific heat 
Cp and isobaric thermal expansion coefficient «p, e.g., see Figs. fl2]) and ([3]). For a water 
monolayer with JV = 1.6 x 10 5 cells confined within nanoparticles with R = 1.6 nm at 
c = 25%, we find a maximum K™ &yi along the isobar at P ~ 0.16 GPa that is 99.7% smaller 
than the c = case. If we decrease c to 2.4%, the reduction of i^™ ax is still remarkable: 
92.3% (Fig. E]). We find similar reductions for C™ ax and a£ ax . 

Such a dramatic reduction of i^™ ax at low T and high P suggests a possible change in 
the region of the phase diagram where water at c = displays the LLPT. From the general 
theory of finite size scaling, we know that at a first-order phase transition i^™ ax , Cp ax and 
Qp as increase linearly with the number of degrees of freedom, here equal to 47V. We find a 
linear increase for 0.14 GPa< P < 0.20 GPa at c = 0, and only for 0.14 GPa< P < 0.16 GPa 
at c = 25% and 2.4%, consistent with the absence of a first-order LLPT outside these ranges. 

To better understand this new feature, i.e., the effect of confinement on the LLPT at high 
P, we study the finite size scaling of the Binder cumulant {23]] U jf = 1 — [(V A ) iy y /3(V 2 } 2 r ], 
where (-)^v stands for the thermodynamic average for a system with jV cells. For jV — > oo, 
at fixed c and P, Ujy = 2/3 for any T away from a first-order phase transition, while 
U™ n < 2/3 at a first-order phase transition 23]. 

For c = 0, we find that Uj n < 2/3 for JY oo at 0.14 GPa< P < 0.20 GPa, while 
U™ n = 2/3, within the error bar, at P = 0.12 GPa (Fig. Hk.). Hence, this analysis confirms 
that for c = there is a first-order LLPT in the range 0.14 GPa< P < 0.2 GPa. 
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For c = 2.4% and 25%, we find that, for large , Uj n < 2/3 at 0.14 GPa, but not 
at 0.12 GPa or at P > 0.16 GPa (Fig. Hb,c). Hence, for c = 2.4% and 25% the first-order 
LLPT occurs only in a limited range of pressures around 0.14 GPa, consistent with our 
results for (5V 2 ) (Fig. [2]) or K™ ax (Fig. [3]), with two end-points: one at ~ 0.15 GPa, another 
at ps 0.13 GPa (Fig. ED . 

We interpret our findings as follows. As a consequence of the stronger HB in the hydration 
shell of each solute, at low T the hydration water is more ordered with respect to the c = 
case. However, shells around different nanoparticles have a different local orientational 
order. This generates competing domains, reminiscent of the locally structured regions 
proposed in Ref. 15J, and exhibits no macroscopic order (upper inset in Fig. [T]). The large 
decrease in fluctuations and response functions, such as Kt, is due to the presence of many 
domain boundaries. Our results for c as low as 2.4% indicates that the decrease is due to 
the introduction of a characteristic length scale, inversely proportional to c, that limits the 
growth of the ordered structured regions. This is consistent also with the results for K^ ax 
(Fig. [3]), where the lower is c, the larger is N beyond which the confined behavior deviates 
from the c = case. 

In previous theoretical analyses, with water confined by a fixed matrix of randomly 
distributed Lennard- Jones disks, the reduction of compressibility was observed only for 



ll| . Here, instead, we find that is reduced 



large hydrophobic obstacle concentrations 
for very low c. 

Our results are qualitatively consistent with recent experiments on H2O confined in the 
hydrophobic mesoporous material CMK-1-14 consisting of micrometer-sized grains, each 
with a 3-dimensional interconnected bicontinuous pore structure, with an average pore di- 
ameter 14 A, at a hydration level of 99% at ambient pressure If]]. Zhang et al. find that the 
TMD is shifted down by 17 K with respect to the hydrophilic confinement in silica mesopores 
and that ap shows a much broader peak, spanning from 240 to 180 K, in contrast to the 
sharp peak at 230 K in hydrophilic confinement 19J], reminiscent of our results on the shift 
of TMD and the reduction of the response functions with respect to the c = case. 

Recent results for small angle x-ray scattering for aqueous solutions of amphiphilic 
tetraalkyl-ammonium cations at ambient conditions suggest that the strengthening of the 
structure of hydration water is present only for solutes with radius smaller than ps 0.44 nm 



241 ] . We therefore repeat our analysis for small nanoparticles with R = 0.4 nm, and find 
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that our results are robust if the amount of hydrophobic interface in contact with water is 
kept constant with respect to the case of R = 1.6 nm. 

In conclusion, we predict that a water monolayer confined in a fixed matrix of hydrophobic 
nanoparticles at concentration c displays changes in the thermodynamics and a drastic 
reduction, > 90%, in Kt, Cp, and otp with respect to the c = case. At c as small as 2.4% 
the first-order LLPT at high P is no longer detected. 
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support by NSF grants CHE0908218 and CHE0911389. GF thanks the Spanish MICINN 
grant FIS2009-10210 (co-financed FEDER). 
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FIG. 1: P-T phase diagram for different nanoparticle concentrations c. Open circles estimate 
liquid-to-gas spinodal line, squares estimate TMD line. In this and all other figures, where not 
shown, errors are smaller than the symbol size. Lines are guides for the eyes (dashed for c = 0, 
dotted for 2.4%, full for 25%). Critical points are shown as large full circles. The liquid-gas critical 
point is the same for c = and 2.4%, while occurs at lower P and T for c = 25%. Lower inset: 
enlarged view of the low-T region. The first-order LLPT ends in a critical point at T ~ 174 K and 
P ~ 0.13 GPa for all c. At c = 2.4% and 25% at P > 0.15 GPa the first-order LLPT is no longer 
detected, indicating a new high-P critical point. Upper inset: configuration at T ~ 160 K and 
P = 0.18 GPa for c = 25%. Hydrophobic nanoparticles are in white; HBs are in different colors 
for different ordered domains. 



8 




Pressure [GPa] 

FIG. 2: (a) Volume fluctuations {5V 2 ) for c = 25% and jY = 10 4 have maxima that follow a locus 
in the P - T plane that does not change, within the error bars, with c or jV . The projections 
{SV 2 ) vs P or vs T clarify that the maxima do not change monotonically with P or T. (b) The 
projection of maxima of (SV 2 ) increase approaching P = 0.132 GPa and 0.156 GPa, consistent 
with our estimate of two critical points at ~ 0.13 GPa and ~ 0.15 GPa. Dashed lines are guides 
for the eyes. 
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FIG. 3: Maxima K^ ax of the isothermal compressibility Kt = (5V 2 ) / {ksT{V)) vs number of water 
molecules N for c = 0, 2.4% and 25%. (a) Linear increase in ET™ ax with N for P = 0.14 GPa, 
consistent with a first-order LLPT for all c. (b) At P = 0.16 GPa, iT™ ax increases linearly for c = 
indicating a first-order LLPT, but saturates for c = 2.4% and 25%, consistent with the absence of 
a first-order LLPT |2J]. 
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FIG. 4: (a) At c = 0, for jV oo is U^ m = 2/3, within the error bar, for P = 0.12 GPa and 
tends to a value < 2/3 for P > 0.14 GPa, indicating a first-order LLPT for P > 0.14 GPa. At 
nanoparticle concentrations c = 2.4% (b) and 25% (c), for .jV — > oo we find ?7™ m < 2/3 only for 
P = 0.14 GPa, indicating that the first-order LLPT is washed out by the hydrophobic confinement 
at high P. For sake of clarity, typical error bars are indicated only for a few points. Lines through 
the points are fits, while other lines are linear interpolations between fits at intermediate P. Black 
arrows mark isobars crossing the first-order LLPT line. 
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